function [R1,R2,R]=tsk_type10(X,M,sigma,C);
R1=[];
R2=[];
[L,n]=size(X);
[m,n]=size(M);

for i=1:L
U=[];
for j=1:m
u=1;
for t=1:n
u=u*gaussmf(X(i,t),[sigma(j,t),M(j,t)]);
end
U=[U,u];
end

c0=C(:,1);
for t=1:n
c0=c0+C(:,t+1)*X(i,t);
end

s=zeros(m,1);

h=U;
delta=U-h;
[l_out,r_out] = interval_wtdavg(c0',s',h,delta);
R1=[R1,l_out];
R2=[R2,r_out];
end

R=(R1+R2)/2;
